Swimming Efficiency of Bacterium Escherichia Coli 
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We use measurements of swimming bacteria in an optical trap to determine fundamental properties 
of bacterial propulsion. In particular, we determine the propulsion matrix, which relates the angular 
velocity of the flagellum to the torques and forces propelling the bacterium. From the propulsion 
matrix, dynamical properties such as forces, torques, swimming speed and power can be obtained 
by measuring the angular velocity of the motor. We find significant heterogeneities among different 
individuals even though all bacteria started from a single colony. The propulsive efficiency, defined 
as the ratio of the propulsive power output to the rotary power input provided by the motors, is 
found to be 0.2%. 



Bacteria swim by rotating helical propellers called flag- 
ella. In the case of Escherichia coli (E. coli), each flagel- 
lum is several microns in length, 20 nm in diameter and 
four to five of them organize into a bundle. The flagella 
are driven at their bases by reversible rotary engines that 
turn at a frequency of approximately 100 Hz. Existing 
experiments show that these molecular engines are Pois- 
son stepping motors consisting of several hundred steps 
per revolution 0, However many essential properties 
of bacterial swimming have not been measured, partic- 
ularly in intact cells. For example, what is the relation 
between the angular velocity of the propellers and the 
force propelling the bacteria forward? What fraction of 
the flagellar motor power is converted into translational 
motion? What variability is there in the swimming ap- 
paratus from cell to cell? Some of these fundamental 
questions have been addressed in theoretical and numer- 
ical work [la- LJ, |5( . however, direct measurements of in- 
tact cells with functional motors and flagella are limited 




Herein, we report an investigation of the fundamen- 
tal swimming properties of E. coli using optical tweezers 
and an imposed external flow. We measure the force re- 
quired to hold the bacterium, and the angular velocities 
of the flagellar bundle and of the cell body as a func- 
tion of the flow velocity. The propulsion matrix, which 
relates the translational and angular velocity of the flag- 
ella to the forces and torques propelling the bacterium, 
can thus be determined one bacterium at a time. Our 
measurements show that although the population aver- 
aged matrix elements are in reasonably good agreement 
with the resistive force theory for helical propellers Q, 
there is a large variability among bacteria from a single 
colony. The propulsion matrix also allows us to deter- 
mine the propulsive efficiency e, defined as the ratio of 
the propulsive power output to the rotary power input, 
to be 0.2%. This is consistent with experiments on he- 
lical propellers Q and close to the maximum efficiency 
for the given cell-body and shape of the flagella. It is 
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FIG. 1: Two different trapping configurations are possible, 
(a) The bacterium can be trapped at the tail of the cell body 
in the presence of an imposed flow U . The trapping is stable 
for U > 0. (b) The bacterium can also be trapped at the head 
of the cell body for U between —40 fim/s and —100 fim/s. 
The forces and velocities are positive if they are along -\-Z. 
The rotations are defined by the right-hand rule such that 
uj < and fl > as depicted, (c) A schematic of a helical 
flagellum: I is the length, 2r is the diameter of the filament, 
9 is the pitch angle of the helix relative to the swimming axis, 
and A is the pitch. 



smaller than the 1 or 2% predicted theoretically for sim- 
ple shapes such as a corkscrew Q. Our experimental 
technique is versatile and can be used to make compar- 
ative studies of mutants strains of the same species or of 
different micro-organisms. Such measurements can shed 
new light on how this remarkable ability to swim evolves 
among different bacterial species. 

An important feature of bacterial swimming is that 
at very low Reynolds numbers (Re ~ 10 -4 ), the fluid 
motion is governed by Stokes flow and nonlinearities in 
the full hydrodynamic equation are irrelevant. Despite 
this simplifying feature, the problem remains theoreti- 
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cally difficult due to complicated time-dependent bound- 
ary conditions. Theoretical studies, therefore, usually as- 
sume that the flagella have very simple geometries such 
as an infinite sheet 0] or a helical coil 0, @ • A second 
approach is not to take into account specific geometries 
but to consider general relations appropriate in the low 
Reynolds number limit 0- In this regime, the torque 
Nfi acting on the propeller and the thrust force Fthrust 
generated by the propeller are linearly related to the pro- 
peller's angular velocity uj and the translational velocity 
v (relative to the background fluid): 



Fthrust = Av-Bu 
N f i = -Bv + Du. 



(la) 
(lb) 



The above equation can be expressed in terms of the sym- 
metric propulsion matrix P ' ^ 



,, ,, , also known 

as the resistance matrix [lfj. Choosing the coordinate 
system in Fig. ^ Fthrust and v are positive if directed 
toward the head of the cell while the sign of u> and Nn 
obeys the right-hand rule, i.e., the flagella is a left-handed 
helix. Based on this coordinate system, the coefficients 
A, B, and D are positive, proportional to fluid viscosity 
?y, and depend on the shape and size of the propeller. The 
basic physics is that in the absence of an applied torque, 
a translating propeller under the influence of an external 
force must rotate, and in the absence of an applied force, 
a rotating propeller under the influence of an external 
torque must translate The above formulation is ap- 
plicable to propellers of any shape and size. However, for 
the special case of a helical coil, the matrix elements can 
be derived from resistive force theory p|: 
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where £ is the length of the coil, and (3 = cos 2 ^> with VP 
being the pitch angle relative to the swimming axis (see 
Fig.nji). The quantity 7 fc is the ratio of the tangential vis- 
cous coefficient K t = 4nr]/(2 ln(0.18A/r) — 1) to the per- 
pendicular viscous coefficient K n = 8^77/(2 ln(0.18A/r) + 
1), where A is the pitch and r is the radius of the coil 
filament. For a smooth coil, Lighthill [4| predicts that 
7fc = K t /K n w 0.7. As can be seen, the helix loses 
its ability to propel, if 7/. — * 1, ^ — * (/3 — * 1) or 
^ — > ? (A — * 0) as expected. 

To complete the model of the swimming bacterium, we 
need the propulsion matrix P for the cell body. Unlike 
P for the flagellum, Pq is diagonal (Bo = 0) since the 
cell body cannot propel itself. The non-viscous force on 
the cell body consists of two parts, the trapping force 
Ftrap due to the optical tweezer holding the bacteria and 
the thrust Fthrust generated by the flagella. The sum of 



these forces should balance the viscous force A Q v acting 
on the cell body. Likewise, the non-viscous torque acting 
on the cell body —Nfi should be balanced by the viscous 
rotational drag. This gives: 



F 



trap 



■ Fthrust — A V, 

D n = -N f i, 



(3a) 
(3b) 



where is the angular velocity of the cell body. We 
treat the cell body as a prolate with minor semi-axis a 
and major semi-axis b so that the linear and rotational 
drag coefficients are A — Anr] fe/(ln(^) — |) and D — 
167i77a 2 i)/3 The optical trapping force is harmonic 

Ft ra p(z) = —k(z — zo), where k is the spring constant 
and z — zq is the displacement from the center of the trap 
[L?! Il3l | . Since the bacteria held by the optical tweezer, 
its net velocity in the lab frame is zero (v' = v + U 0), 
and the relative velocity v is opposite to the external flow 
U. Substituting v = —U into Eqs. Ill and l3l gives: 

k{z-z ) = [A + A )U + Buj, (4a) 
D Q = -BU-Duj. (4b) 

This set of equations will be used below to analyze our 
data. 

We used a non-tumbling strain of bacteria (RP5231) 
in our measurements. We were delighted to find that 
such a bacterium swimming a few microns above a glass 
surface could be stably trapped, along its swimming di- 
rection, by the optical tweezer |14| . The bacterium can 
then be manipulated by an imposed external flow. Fig- 
ure n illustrates our experimental setup along with the 
flow configurations. A bacterium swimming to the left 
(along the +Z direction) is held by a strongly focused 
IR laser (X— 1024 nm). In the absence of the flow, the 
bacterium is invariantly held by the tail of the body. The 
thrust force and the trapping force are balanced and the 
bacterium is stationary with respect to the trap. The 
trapping remains stable when a uniform flow in the +Z 
direction is applied (U > 0) [l^. The bacterium can also 
be trapped at the head of the body, but the flow field 
must be reversed (U < 0). 

To measure the trapping force and the position of the 
trapped cell tip, the transmitted IR beam was refocused 
by a high numerical aperture condenser (N.A. 1.5) and 
projected onto a two-dimensional position sensitive de- 
tector (Pacific Silicon Sensor Inc., DL100-7PCBA). The 
position of the trapped cell tip with respect to the cen- 
ter of the trap is monitored by a PC equipped with a 
National Instruments analog-digital converter card (AT- 
MIO-16E-2). The conversion rate is 10 kHz at 12-bit 
resolution. A non-flagellated bacterium was used to cal- 
ibrate the spring constant k of the optical trap by mea- 
suring the position of the trapped tip as a function of 
the flow U. For the laser intensity (23 mW) used in this 
experiment, k = 5.7 x 10~ 6 Af/m. A brief description 
of the calibration process is presented in Materials and 
Methods. 

Figure|21displays an example of the time trace z(t), the 
longitudinal displacement along the swimming direction 
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FIG. 2: A typical experimental run for a swimming bacterium 
held in the optical trap. An uniform flow U = 10 fim/s 
is established in regime I. The flow U is decreased to zero 
linearly in regime II. The flow U remains zero in regime III. 
The laser is switched off momentarily to let the bacterium 
escape and the undeflected laser beam position is recorded in 
regime IV. Solid lines depict linear fits to each regimes. 



of a trapped bacterium. We observed large oscillations 
overlying a systematic variation of z(t) as the external 
flow is changed. These oscillations result from wobbling 
of the cell body in response to the rotation of the flagella 
bundle [Til IT?! ] . The trapped bacterium was perturbed 
by the following sequence of events: In regime I, the bac- 
terium is subject to a uniform flow U = +10 fim/s. The 
bacterium maintains an average position away from the 
center of the trap. In regime II, U is linearly reduced to 
zero in 5 s. The average bacterium position shifts system- 
atically toward the center of the trap. In regime III, U is 
maintained at zero for 5 s, and the average position of the 
bacterium relative to the trap is again constant. Finally 
in regime IV, the bacterium is released. The position 
of the undeflected beam in regime IV is taken to be z , 
the center of the optical trap. From regime II we obtain 
the net translational drag coefficient A + Aq = kAz/AU, 
and in regime III we obtain F thrust , since F trap = F thrust 
when [7 = 0. We checked that this measurement was re- 
producible by returning the flow to U = 10 /im/s rather 
than releasing the bacterium after regime III. The bac- 
terium returned to within a few percent of its initial av- 
erage position. 

We extracted the angular velocities and u> using a 
Fourier analysis of the time trace x(t) of the transverse 
position of the cell body. This transverse signal shows 
more pronounced oscillations than z{t). Figure E^a) dis- 
plays a sample power spectrum E(f) for a short time 
trace of 5 s when [7 = (regime III). The power spectrum 
has two strong peaks at /l — 7 Hz and Jh — 130 Hz, 
respectively. These two frequencies can be associated 
with the angular velocities of the cell-body SI = 2irfL 
and of the flagellum bundle oj = — 27r/# 16]. Averag- 
ing over 250 bacteria, we found Jl = (8.0 ± 0.2) Hz and 
Jh = (125 ± 2) Hz, where the ± are standard errors 
of the mean. However, as shown in Fig. 0tc-d) there 
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FIG. 3: (a) Power spectrum of x(t) shows peaks correspond- 
ing to /l and fn. (b) The variation of the rotation frequency 
of the cell body fa as a function of flow speed U. The linear 
dependence is consistent with the propulsion matrix formu- 
lation, (c) and (d) delineate the PDFs of fa and fn, re- 
spectively. Error bars are standard errors of the mean unless 
otherwise noted. 



is considerable variation of and fjj between individ- 
ual bacteria. The standard deviations er/ L = 2.4 Hz and 
<jf H =27 Hz are respectively 20 and 30 % in the mean 
values. 

To test the basic physics implied by the propulsion ma- 
trix, we measured the dependence of fa and fjj on U for 
an additional 40 bacteria which were subjected to flow 
speeds of [7 = —40,-60,-80 and — 100_/xm/s. Figure 
|3Jb) shows that the average frequency fa increases lin- 
early with —[7 and the result is in good agreement with 
Eq. I4bl as predicted by the propulsion matrix formula- 
tion. Within the noise of the measurement, no system- 
atic change in fu was detected. This is expected since 
at low-loads the molecular motor is known to rotate at a 
constant angular speed independent of the load 0] . 

To complete our determination of the propulsion ma- 
trix, the semi-minor axis a and length L ce u = 2b of 
the bacteria were measured directly by video microscopy 
while immobilized in the trap. This allows us to calcu- 
late the drag coefficients Ao and Dq for the cell body. 
From the time trace z(t), A and B are calculated by 
A = kAz/AU - A and B = F thrust /w when [7 = 0. 
Finally, the measurements of the angular velocities gives 
D = —Dqlu/Q. The matrix elements averaged over a pop- 
ulation of 250 bacteria are A = (3.8 ±0.2) x 10~ 8 N s/m, 
B = (4.0 ± 0.1) x 10~ 16 TV s, and D = (5.9 ± 0.1) x 
10~ 22 7V s m. The translational drag coefficient of the 
flagella is approximately twice that of the cell body 
(Ao = 1.7 x 10~ 8 7V s/m). Therefore a significant portion 
of drag is due to the flagella. On the other hand, the ro- 
tational drag of the flagella D is much smaller than that 
of the cell body (D = 8.8 x 10" 21 7V s m). 

It is instructive to use the measured propulsion matrix 
to extract physical parameters that are relevant to flag- 
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FIG. 4: The PDFs of (a) A/A, (b) B/B, and (c) D/D. The insets shows the length L ceU dependence of (a) (A), (b) (B), and 
(c) (D). The solid lines in the main figures are fits to log- normal distributions. 



ellar bundles. The resistive force theory for the helix coil 
contains four independent parameters: the pitch A, the 
pitch angle 'F, the length of the helix I, and the radius r 
of the filament, assuming that the viscosity of the fluid 
(r/ = 10~ 3 Pa s) is known. The three matrix elements A, 
B, and D in Eqs. [2] however are not sufficient to predict 
all the four geometric parameters. To make progress, we 
assumed that the the pitch angle 'F is 41°, as determined 
by Turner et al. using fluorescently labeled E. coli cells 
0- This angle also turns out to be remarkably close to 
the optimal angle (42°) that maximizes the propulsion ef- 
ficiency of an ideal helix Using /3 = cos 2 (41°) = 0.57, 
Eqs. [2] predict = 0.84, A = 0.9 ym, r — 23 nm, and 
£ = 6.2 ym. These values are comparable to the fluores- 
cent measurements of Turner et al. who found £ = 7 ym, 
and A ~ 1 ym for the curly flagella and 2.2 ym for the 
normal ones 7] . We can estimate the average number of 



flagella using r 



<N 



which gives N 



A 



5.3. 



This is slightly greater than the N sa 3.3 found by Turner 
et al. The difference may be expected because the bac- 
teria used in their experiment are shorter than the ones 
we studied; longer bacteria usually have more flagella. 

All important dynamical quantities can be obtained 
from our measurements. For example, the average thrust 
for U — is Fthrust — B Hi = 0.31 pN. The average 
swimming speed is V SW i m = B u>/(A Q + A) = 6 ym/s, 
which should be compared with V SW im ~ 10 ym/s we ob- 
tained directly by video microscopy. The difference may 
be due to correlations between A, B and D; both B and 
D grows with A on average. Similarly the average torque 
Nfi — D (\Q — Q,\) = 4.9 x 10~ 19 N m is surprisingly close 
to that found for Streptococcus |l8j . 

We observed that the propulsion matrix elements vary 
greatly among individual bacteria even though our bac- 
teria are from a single colony. Figure 0] shows the prob- 
ability distribution functions (PDF) of the scaled quan- 
tities A/A, B/B, and D/D. The standard deviations 
a are comparable to the mean values with o~a/A = 0.7, 
ctb/B = 0.5 and oj^/D = 0.3. A conspicuous feature of 
all the PDFs is their broad tails, particularly for the lin- 
ear drag coefficient A. This might be an indication of ei- 
ther significant structural heterogeneities in the flagellar 
bundles of individual cells or of changes in conformations 
of the flagellar bundle with time. As is often the case in 



biological systems, the PDFs can be roughly fitted to log- 
normal distribution functions, which are plotted as solid 
lines in the figure. 

Part of the variation in A, B and D must arise due 
to bacteria being in different stages of their growth cy- 
cle during the measurements. This is particularly the 
case in the middle log-phase of a growing culture, where 
the bacterial size is highly varied. For ease of trapping, 
very long and very short bacteria were excluded from the 
measurements and the middle-sized bacteria (4—12 ym), 
which comprised about 95% of the population, was cho- 
sen. The bacterial cell- length distribution of this selected 
population is plotted in Fig. [Sfa). The figure also shows 
the cell length at which we first observed septal rings 
(dotted line) and the length at which cells divide (solid 
line). We used the bacterial length L ce u as a measure 
of the bacterium's physiological state and determine the 
propulsive matrix elements as a function of L ce u . To de- 
termine the length dependence of the coefficients A, B, 
and D, we calculated the averaged values (A), (B) and 
(D) for bacteria of similar length. The result are pre- 
sented in the insets of Fig. 0Ja-c). The linear drag coeffi- 
cient (A) has no clear size dependence but (B) is peaked 
at L ce u sa 6 ym, which coincides with the peak of the 
size distribution. On the other hand, Fig. 0Jc) shows 
that the rotational drag coefficient (D) of the propeller 
grows linearly with the cell-body length L ce u. 

These size dependencies allow us to draw conclusions 
about the structure of flagellar bundles at different stages 
of cell growth. Inspection of Eqs. [5] shows that the ma- 
trix elements scale with the pitch A according to A oc A , 
B « A 1 , and D oc A 2 . The fact that we find that B 
and D depend on L ce u but A does not, implies that the 
primary L ce u dependence is in pitch A. The pitch angle 
^ and the flagella length £ are approximately constant 
independent of L ce u. Likewise, 7^ depends logarithmi- 
cally on A and so is only very weakly dependent on L ce u . 
Since our measurements show a linear relation between 
D and L ce u , it suggests that A 2 grows linearly with L ce u . 
A possible scenario for this is that more and more flag- 
ella are incorporated into the bundle as the bacteria cell 
body grows and this causes the A 2 to grow in agreement 
with Fig.^Jc). From the shortest to the longest bacterial 
body size (L ce u), we found that the fractional change 
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6X/X is about 16%, which may be discernible in care- 
fully conducted observations using fluorescently labeled 
bacteria. 

We next turn our attention to the power and propulsive 
efficiency of the swimming bacteria. The average power 
output of the flagellar motors is E = Dofl(oj — fi) = 
0.4 pW. The power used to turn the cell body is 
Dofl 2 w 0.02 pW while the actual propulsive power is an- 
other factor of ten smaller with A V 2 wim w 0.0017 pW. 
Therefore 5% of the rotary power is used to rotate the 

I 



Here we used B 2 < (A a + A)D and D > D to obtain 
the approximate form. These assumptions are met on 
average but does not always hold for a particular bac- 
terium. Therefore we use the full form to calculate the 
efficiency. FigureOJc) shows that the efficiency as a func- 
tion of bacterial size is constant up to the cell division 
length (L ce ii = 8 p,m). The average efficiency e » 0.2% 
(l9| is surprisingly close to that of sedimenting helices 
in a silicon oil, which were tested as model flagella by 
Purcell |3 ■ We can also ask, for a given A , what is 
the maximum efficiency attainable by the bacterium as 
a function of the length of the flagella. Assume that at 
some characteristic length £ p , the propulsive coefficients 
of the flagellum are A p , B p and D p . Neglecting logarith- 
mic corrections and assuming the width of the flagellar 
bundle is constant, these coefficients should grow linearly 
with flagella length £ so that A w kA p , B w kB p and 
D « kD p where k = £/£ p . This assumption is consis- 
tent with Eqs. |2j Substituting for A, B and D into 
our approximate expression for e (in Eq. [5J wc find that 
the maximum efficiency occurs when A = Aq and that 
£max ~ B p /(4:A p D p ) which depends only on the shape of 
the propeller. The same result was obtained by Purcell 
when he maximized e by assuming that all propeller di- 
mensions (not just the length) scaled with n [9j. In our 
experiments, we find that A w 2Ao and so flagella are 
twice as long as that required to maximize its propul- 
sive efficiency. However, the peak in e as a function of n 
is fairly broad and the observed efficiency is about 75% 
of the maximum efficiency as shown by the solid and 
dashed lines in Fig. OJ. The broadness of the peak may 
be why the propulsive efficiency is approximately con- 
stant throughout the bacterial cell division cycle. The 
experimental propulsive efficiencies are consistent with 
the small e max of between 0.3% — 0.8% measured for he- 
lical propellers 0- It is smaller still than the 1 or 2% 
predicted theoretically for a helical propeller 0, . 

In summary, bacterial propulsion in a uniform stream 
is investigated with the help of optical tweezers, which 
allow the thrust force Fthrust to be directly measured as 



cell body, and only 0.5% is used to push the bacteria for- 
ward. Figure EJb) shows the average motor power as a 
function of bacterial length L ce u. The power increases 
gradually with L ce u, which is consistent with the above 
discussion that the number of flagella N and the associ- 
ated motors increase with L ce u . The propulsion efficiency 
e, defined as the ratio of the propulsive output power to 
the rotary input power, can be related to the propulsion 
matrix elements 91: 



(5) 



a function of imposed flow. For a free swimming bac- 
terium, Fthrust precisely balances the viscous drag of the 
cell body A$v and of the flagellar bundle Av. The contri- 
bution to the drag by Av is twice as large as A a v but is 
difficult to determine without our direct force measure- 
ments. We also showed that the propulsion matrix de- 
scription proposed by Purcell gives an adequate descrip- 
tion of bacterial propulsion over a physiological range of 
velocities. In retrospect, the validity of the propulsion 
matrix, or for that matter the resistive force theory it- 
self, is not self evident for real micro-organisms because 
of possible deformations of flagellar bundles due to hy- 
drodynamic stresses induced by swimming or by the flow 
5J- 

We have determined all elements of the propulsion ma- 
trix and used the resistive force calculations for a helical 
coil to estimate microscopic properties of flagella . The 
results were consistent with earlier measurements even 
though the resistive force calculations neglect the effect 
of long-range hydrodynamics interactions between differ- 
ent parts of the flagella and with the cell body. Using 
the propulsion matrix, we have determined all dynamic 
quantities related to bacterial swimming and their de- 
pendence on the size of the cell body. In particular, we 
found that the propulsive efficiency e, defined as the ratio 
of the propulsive power output to the rotary power input, 
is ~ 0.2% and is nearly independent of L ce u. The mea- 
sured e is consistent with the experiments on rigid helical 
propellers [t| and close to the maximum efficiency for the 
given size of the cell body and the shape of the flagel- 
lar bundle. A conspicuous finding of our measurements 
is that all the matrix elements are broadly distributed 
despite the fact that all bacteria started from a single 
colony. The ubiquity of such broad distributions in bio- 
logical systems is significant and begs further systematic 
study. 



MATERIALS AND METHODS 

Sample preparation: We followed standard growth con- 
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FIG. 5: (a) The PDF of the bacterial cell length L ce u- The dashed line is a fit to the log-normal distribution. The dotted and 
solid vertical lines are, respectively, the cell lengths at which we first observed a septal ring and where cell division occurred, 
(b) The flagellar power output (E) as a function of L ce u. The dashed line is a linear fit. (c) The propulsion efficiency (e) as a 
function of L ce u. The dotted horizontal line marks the mean efficiency 0.2% of the entire population. The solid horizontal line 
is the maximum efficiency e max when the flagellum length is optimized. See text for details. 



ditions for culturing bacteria E. coli strains, RP5231 
and YK4516. RP5231 is a smooth swimming strain be- 
cause two of its chemotactic genes, CheY and CheZ, 
were deleted. A single colony was picked from a fresh 
agar plate and grown to saturation overnight in the LB 
medium (peptone 4 g, yeast 2 g, NaCl 1 g, 1M NaOH 0.4 
ml; for 400 ml of media). The culture was maintained 
at 30 °C and was shaken continuously at 200 rpm. The 
overnight sample was diluted 1:100 in fresh LB medium 
and grown to the middle log phase for 3 hours. 

To calibrate the spring constant k of the optical 
trap, we used non-flagellated bacteria strain YK4516. 
A uniform flow U was applied and the shift in the 



centroid of the transmitted IR laser beam was recorded 
by the position detector. For an ellipsoid body, the 
translational drag coefficient A is known and the spring 
constant is obtained using k — A AU/ Az. The noise 
in the output of the optical trap was 0.1 nm/^/Hz (for 
and sampling rate lOKHz). 
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